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Nonlinear wave dynamics near phase transition in TT-symmetric localized potentials 
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Nonlinear wave propagation in parity-time {VT) symmetric localized potentials is investigated 
analytically near a phase-transition point where a pair of real eigenvalues of the potential coalesce 
and bifurcate into the complex plane. Necessary conditions for phase transition to occur are derived 
based on a generalization of the Krein signature. Using multi-scale perturbation analysis, a reduced 
nonlinear ODE model is derived for the amplitude of localized solutions near phase transition. 

Above phase transition, this ODE model predicts a family of stable solitons not bifurcating from 
linear (infinitesimal) modes under a certain sign of nonlinearity. In addition, it predicts periodically- 
oscillating nonlinear modes away from solitons. Under the opposite sign of nonlinearity, it predicts 
unbounded growth of solutions. Below phase transition, solution dynamics is predicted as well. All 
analytical results are compared to direct computations of the full system and good agreement is 
observed. 


I. INTRODUCTION 

Parity-time {VT) symmetric systems started out from 
an observation in non-Hermitian quantum mechanics, 
where a complex but TT-symmetric potential could pos¬ 
sess all-real spectrum @|. This concept later spread out 
to optics, Bose-Einstein condensation, mechanical sys¬ 
tems, electric circuits and many other fields, where a 
judicious balancing of gain and loss constitutes a VT- 
symmetric system which can admit all-real linear spec¬ 
trum I2j4l2j. For example, in optics an even refractive 
index profile together with an odd gain-loss landscape 
yields a TT-symmetric system. A common phenomenon 
in linear PT-symmetric systems is phase transition (also 
known as TT-symmetry breaking), where pairs of real 
eigenvalues collide and then bifurcate to the complex 
plane when the magnitude of gain and loss is above a 
certain threshold 0,1 EMI- This phase transition has 
been observed experimentally in a wide range of physical 
systems @1,0, Ho . 11 1 . When nonlinearity is introduced 
into VT systems, the interplay between nonlinearity and 
VT symmetry gives rise to additional novel properties 
such as the existence of continuous families of station¬ 
ary nonlinear modes, stabilization of nonlinear modes 
above phase transition, and symmetry breaking of non¬ 
linear modes [(I |l~4l422l . These findings reveal that PT- 
symmetric systems break the boundaries between tra¬ 
ditional conservative and dissipative systems and open 
new exciting research territories. Practical applications 
of VT systems are starting to emerge as well, such as re¬ 
cent demonstrations of PT-symmetric micro-ring lasers 
and unidirectional VT metamaterials @, Hi, [J3, Ed). 

An important feature of PT-symmetric systems is 
phase transition, where the linear spectrum changes 
from all-real to partially-complex and infinitesimal waves 
change from stable to unstable. At phase transition, a 
pair of real eigenvalues coalesce and form an exceptional 
point with a non-diagonal Jordan block (i.e., with the 
algebraic multiplicity higher than the geometric multi¬ 
plicity). Phase transition is a distinct linear property of 
PT-symmetric systems, and it is at the heart of many 
proposed applications such as PT-symmetric micro-ring 


lasers and unidirectional VT metamaterials j§j, IHJ . 

When nonlinearity is present (such as if the wave am¬ 
plitude is not small), the interplay between phase transi¬ 
tion and nonlinearity is a fascinating subject. This inter¬ 
play was previously studied for periodic PT-symmetric 
potentials in [21], j2J, (25|, where novel behaviors such as 
wave-blowup and oscillating bound states were reported 
below phase transition. In addition, stable nonlinear 
Bloch modes were reported above phase transition be¬ 
cause nonlinearity transforms the effective potential from 
above to below phase transition [21,]. However, in peri¬ 
odic potentials above phase transition, the presence of 
unstable infinitely extended linear modes makes the zero 
background unstable, which excludes the possibility of 
stable spatially-localized coherent structures. In local¬ 
ized potentials, will the situation be different? 

In this article we study nonlinear wave behaviors in lo¬ 
calized PT-symmetric potentials near phase transition. 
Unlike periodic potentials, the instability of linear modes 
above phase transition is limited to the area around the 
localized potential. In this case, the addition of nonlin¬ 
earity can balance against gain and loss making stable 
spatially-localized coherent structures, such as solitons 
and oscillating bound states, possible above phase tran¬ 
sition. Mathematically we explain this phenomenon by a 
multi-scale perturbation analysis, where a reduced non¬ 
linear ODE model is derived for the amplitude of local¬ 
ized solutions near phase transition. Above phase transi¬ 
tion, this ODE model predicts a family of stable solitons 
not bifurcating from linear (infinitesimal) modes under a 
certain sign of nonlinearity. In addition, it predicts per¬ 
sistent oscillating nonlinear modes away from solitons. 
Under the opposite sign of nonlinearity, it predicts un¬ 
bounded growth of solutions. Similarly, solution dynam¬ 
ics below phase transition is predicted as well. All these 
predictions are verified in the full PDE system. In addi¬ 
tion to these nonlinear dynamics, we also derive a neces¬ 
sary condition for phase transition to occur at an excep¬ 
tional point in the linear VT system by a generalization 
of the Krein signature, namely, phase transition from a 
collision of two real eigenvalues is possible only when the 
two eigenvalues have opposite TT-Krein signatures. 
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II. PRELIMINARIES 

The mathematical model we consider in this article is 
the following potential NLS equation 

itl>z + ^xx + V(x-,e)i> + a\tp\ 2 ip = 0, (2.1) 

where V{x\e) is a PP-symmetric complex potential, i.e., 

V*(-x-,e)=V(x-,e), (2.2) 

parameterized by e which controls the gain-loss strength, 
a = ±1 is the sign of nonlinearity, and the superscript 
**’ represents complex conjugation. This model governs 
nonlinear light propagation in an optical medium with 
gain and loss [14J as well as dynamics of Bose-Einstein 
condensates in a double-well potential with atoms in¬ 
jected into one well and removed from the other well Q. 
Without loss of generality, we assume phase transition 
occurs at e = 0, where a pair of real eigenvalues of the 
potential coalesce and form an exceptional point, and we 
will analyze the solution dynamics in Eq. m near this 
exceptional point, i.e., when |e| -C 1. 

The analysis to be developed applies to all localized 
PP-symmetric potentials near phase transition. To il¬ 
lustrate these analytical results and compare them with 
direct numerics of the full model (12.11) , we will use a con¬ 
crete example — the so-called Scarff II potential 

V = VRsech 2 (x) + iWosech(a;)tanh(x), (2.3) 

where Vr, Wq are real parameters. For this potential, 
phase transition occurs at Wq = Vr + 1/4 [13], and soli- 
tons as well as robust oscillating solutions were reported 
numerically below phase transition in [TJ, [pij - lpl . 

III. PP-KREIN SIGNATURE AND A 
NECESSARY CONDITION FOR PHASE 
TRANSITION 

We begin by studying the general linear eigenvalue 
problem 

L(x\ e)u = — fm, (3-1) 

where L is a PP-symmetric linear operator parameter¬ 
ized by e, i.e., 

L*(-x-,e)=L(x-,e), (3.2) 

and /r is an eigenvalue. We wish to consider the phase- 
transition process by which the spectrum of L changes 
from all-real to partially-complex. This phase transi¬ 
tion occurs when a pair of real eigenvalues collide and 
then bifurcate into the complex plane. It is important 
to recognize that not any two real eigenvalues can turn 
complex upon collision. This is analogous to the linear 
stability of equilibria in Hamiltonian systems, where not 
just any two purely imaginary eigenvalues upon collision 


can bifurcate off the imaginary axis and create linear in¬ 
stability f29j - [3lT . Then the question we address is: under 
what conditions can a pair of real eigenvalues of L induce 
phase transition upon collision? 

For the potential NLS equation m, when one looks 
for linear eigenmodes i/> = u(x)e~ llJZ , the eigenvalue 
problem (13.11) will be obtained with 

L = d xx + V(x;e ), (3.3) 

which is PP-symmetric. However, in this section we will 
consider the eigenvalue problem (13.11) for general PP- 
symmetric operators, not just (13.31) . 

First we make some elementary observations. Since 
the operator L is PP-symmetric, any complex eigenval¬ 
ues must come in conjugate pairs, i.e., if [/i,u(x)] is an 
eigenmode, then so is \jj,* ,u*(—x)]. If /r is a simple real 
eigenvalue, then its eigenfunction u(x) can be made PP- 
symmetric by scaling. 

We start the analysis by introducing a sesquilinear PP- 
product 

/ OO 

f*(~x)g(x)dx, ( 3 . 4 ) 

-OO 

which naturally satisfies the symmetry condition 

(gj)rr = (f,g)vr- (3-5) 

Thus for any complex function /( x), (/, f)vr is rea l and 
invariant under a gauge transformation f{x) —> f(x)e ie , 
where 9 is a real constant. In addition, under this PP- 
product the PP-symmetric operator L is “self-adjoint”, 
i.e., 

(Lf,g)vr = (f,Lg)rr- (3.6) 

For an eigenmode [g, u(x)] of (13.31) with a simple real 
eigenvalue g we define its PP-Krein signature as 

S(li) = sgn[(u,u)73 r ]. (3.7) 

Here (u,u)-pT 7^ 0, because u{x) can be made PP- 

symmetric, i.e., u*{—x) = u(x), so ( u,u)'pt = /_ M 
which is nonzero since it is the Fredholm condition for the 
generalized-eigenfunction equation (L + n)u g = u not to 
admit a solution in view that u* ( x ) is in the kernel of the 
adjoint operator L* + g. Thus the PP-Krein signature 
of a simple real eigenvalue is always positive or negative. 
In addition, this signature cannot change under contin¬ 
uous variation of the parameter e unless pairs of such 
eigenvalues collide. 

The main result of this section is that when two such 
eigenvalues collide, a necessary condition for complex- 
eigenvalue bifurcation is the two real eigenvalues have 
opposite PP-Krein signatures. This result extends an 
analogous one in Hamiltonian systems to PP-symmetric 
systems [29l [3lj . 

We first present three lemmas. 
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Lemma 1 Let u\{x) and m 2 (x) be two eigenfunctions of 
the VT-symmetric operator L with real eigenvalues pi 
and p ,2 respectively. If p i 7 ^ p 2 then fui^uffvT = 0. 

Proof Using the “self-adjoint” property (13.61) we have 
(u\,Lu 2 )vt = (Lui, uffvT■ 


quadratic form (/, f)vr restricted to the subspace S = 
span(Mi,M 2 ) by looking at the dual matrix M. 

When e > 0 and pi = p%, we see easily by Lemma [2] 
that 


Me>o — 


0 b 
b* 0 ’ 


(3.9) 


Then using the fact that ui and U 2 are eigenfunctions, 
we can calculate the left and right sides of the above 
equation as 


{u± , Lufd'PP — ~P2(ui,U2)vr> 


and 


(Lui,U 2 )vr = “Ml( u l, u 2 )pT- 

Thus if pi 7 ^ P 2 then {u\,U 2 )vt = 0. □ 

Lemma 2 Let u( x) be an eigenfunction of the VT - 
symmetric operator L with a complex eigenvalue p. Then 
{u,u) VT = 0. 

Proof From Eq. (13.61) we have ( u , Lu)-pj- = (Lu, u)-pr- 
Calculating the two sides of this equation we get 

i.h'- h*){u 1 u) VT = 0 . 

Thus if p is complex, then (u , u)-pT = 0. □ 

Lemma 3 Let {e^} be a basis for an N-dimensional 

N 

functional subspace and f = CjCj. Then 

j=i 

(/) f)vr = c H Mc, 

where M is a N x N Hermitian matrix with elements 
given by 


Mij = (ei,ej)-pTi (3-8) 

and the superscript ‘H ’ represents the Hermitian of a vec¬ 
tor. 

Proof Substituting the / expression into (/, f)vr and 
utilizing the linearity of the PP-product, this lemma can 
be readily proved. The Hermiticity of M comes directly 
from relations (13.51) and arm . □ 

Now we present the main result of this section. 

Theorem 4 Let L be a VT-symmetric operator param¬ 
eterized by e. If a pair of simple real eigenvalues of L 
collide and bifurcate into the complex plane at e = 0 , 
then before the bifurcation the two real eigenvalues must 
have opposite VT-Krein signatures. 

Proof Let ui(x) and u 2 (x) be two eigenfunctions of 
L with eigenvalues pi and P 2 respectively; when e < 0, 
pi and P 2 are simple real with pi ^ P 2 , and when e > 0 , 
pi and P 2 are complex with pi = p%. We analyze the 


where b = (ui,U 2 )vt■ I n this case, ui(x) = u%(—x), thus 

/ OO 

u\{x)dx , 

-OO 


which is nonzero because for a simple eigenvalue P 2 , 
f_ u%dx 7 ^ 0 is the Fredholm condition for non¬ 
existence of a generalized eigenfunction [see earlier text 
below Eq. (EH)]. This means M has a pair of real eigen¬ 
values of opposite sign and is thus indefinite. 

Likewise, for e < 0, since pi and P 2 are strictly real, 
then in view of Lemma [TJ 


M e <o — 


ai 0 

0 a 2 


(3.10) 


where ai = (ui,ui)-pj- and a 2 = ( U 2 ,U 2 )vt have the 
signs of the PT-Krein signatures for pi and /r 2 respec¬ 
tively. 

At e = 0 the two eigenvalues collide, and pi = P 2 = 
po. There are two cases here. The first case is where 
L has two eigenfunctions ui and u 2 at po , i.e., L has 
a diagonal Jordan block. In this case, ( Uk,Uk)vr 7 ^ 0 
(k = 1 , 2 ), and we can always choose the basis {■ui,m 2 } 
so that (ui,U 2 )vt = 0 , thus the operator M e=0 has the 
same structure as (13.101) . The second case is where L has 
a single eigenfunction uq at po, i.e., L has a non-diagonal 
Jordan block. In this case, the subspace S reduces to S = 
span(wo, %), where u g is the generalized eigenfunction 
satisfying 


C Po + L)u g = Mo- 

Taking the PT-product of this equation with uq and re¬ 
calling EH), i-e., applying the Fredholm condition, we 
get (u 0 ,u 0 )vr = 0 , thus 


M e=0 = 


0 ci 

ct c 2 


(3.11) 


where Ci = {uQ,Ug)'pj- and c 2 = (u 9 , M g )-p 7 -. Note that 
Ci 7 ^ 0 since it is the Fredholm condition for eigenvalue 
po not to have a second generalized eigenfunction. Then 
det(M e=0 ) < 0 , hence M e=0 has a pair of real eigenvalues 
of opposite sign and is indefinite. 

Since M is indefinite for e > 0 and this indefiniteness 
is continuous across the bifurcation point e = 0 , we see 
that M for e < 0 must also be indefinite, which directly 
implies that real eigenvalues pi and /i 2 before bifurcation 
must have opposite PT-Krein signatures. □ 

Now we use an example to illustrate this theorem. In 
the Scarff-II potential (12.31) . we fix Vr = 5 and vary the 
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gain-loss coefficient Wq. The linear spectra for various 
W 0 values are displayed in Fig. [T] It is seen that phase 
transition occurs at Wq = 5.25, where a pair of simple 
real eigenvalues coalesce and form an exceptional point, 
which then turns complex when Wo > 5.25. We have 
calculated the PT-Krein signatures of those real eigen¬ 
values (indicated by colors in the figure) and found them 
indeed opposite, in agreement with Theorem [I] 


2 4 6 

Wo 

FIG. 1: Eigenvalues h for varying gain-loss strength Wo in the 
Scarff-II potential (IPll with Vr = 5 (the continuous spectrum 
is displayed in light blue). Blue lines indicate eigenvalues with 
positive PT-Krein signatures while red lines indicate eigenval¬ 
ues with negative PT-Krein signatures. Complex eigenvalues 
are indicated with dotted purple lines. 

Interestingly, Fig. [1] also shows another collision of sim¬ 
ple real eigenvalues of opposite PT-Krein signatures at 
Wo = 4.25 (which creates another exceptional point). 
However, after collision these real eigenvalues re-emerge 
and no complex eigenvalues bifurcate out. This shows 
that collision of real eigenvalues of opposite PT-Krein 
signatures is a necessary but not sufficient condition for 
complex-eigenvalue bifurcation. A sufficient condition for 
complex-eigenvalue bifurcation can be found in later text 
[i.e., «/0, see the paragraph below Eq. (14.121) ]. This 
condition is not met at the other exceptional point, thus 
we do not see phase transition there. 
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IV. REDUCED MODEL NEAR PHASE 
TRANSITION 

In this section, we consider the potential nonlinear 
Schrodinger equation m and analyze its solution dy¬ 
namics near phase transition. 

Let us suppose the PT-symmetric potential in 
Eq. (12.11) takes the form 

V(x;e) = V 0 (x) + e 2 V 2 {x), (4.1) 

where 0 < e ^ 1 . Here Vq(x) is the unperturbed poten¬ 
tial, V 2 (x) is the form of potential perturbation, and e 2 
is the strength of this perturbation. We assume that the 
unperturbed potential Vo(x) is at phase transition and 
possesses an exceptional point at fj, = Ho, be., the linear 
operator d xx + Vq(x) has a single eigenfunction u e (x) and 
a generalized eigenfunction u g (x ) at fj, q . Defining 


(4.2) 


then we have 


L 0 u e = 0, L 0 u g = u e . (4.3) 

Since Lq is PT-symmetric, both u e and u g can be chosen 
to be PT-symmetric as well, i.e., 


u* e {-x) = u e (x), u*(~x) = u g {x). (4.4) 

By taking the complex conjugate of the Lou e = 0 equa¬ 
tion, we see that it* is in the kernel of the adjoint opera¬ 
tor L*. Thus the solvability condition for the Lou g = u e 
equation is that 



(4.5) 


In principle, an exceptional point can have algebraic 
multiplicities higher than two, meaning that it can have 
additional generalized eigenfunctions beside u g . But in 
a generic case, an exceptional point is formed by the col¬ 
lision of two simple real eigenvalues, in which case its 
algebraic multiplicity is only two. For simplicity, we only 
consider such generic exceptional points in this article. 
Since their algebraic multiplicities are two, they do not 
admit other generalized eigenfunctions, i.e., the equation 


Lou g 2 = u g 

admits no localized solutions for u g 2 - Since it* is in the 
kernel of the adjoint operator Lq, the Fredholm condi¬ 
tion on the above equation is that its right side be not 
orthogonal to it*, i.e., 

/ OO 

u e u g dx 0. (4-6) 

-OO 

In addition, since u e and u g are PT-symmetric, so is 
u e u g , hence D is real. 

If the potential Vo ( x ) is perturbed to be dm we study 
nonlinear dynamics in this perturbed potential by multi¬ 
scale perturbation methods below. First we expand the 
solution to Eq. HO into a perturbation series, 

ij){x, z) = (eni(x, Z) + e 2 u 2 + e 3 it 3 + ...) e -1MoZ , (4.7) 

where Z = ez. Then up to order e 3 we have a system of 
equations 


LqU\ — 0, 

LqU 2 = —lUiz, 

Lqu 3 = -iu 2 z - V 2 ui - ct | ui | 2 ui . 

Since it* is in the kernel of the adjoint operator Lq, the 
solvability conditions for these equations are that their 
right sides be orthogonal to the adjoint homogeneous so¬ 
lution it*. 

At orders e and e 2 we find from (14.31) that 

iti = A(Z)u e (x), 
u 2 = -iA z Ug. 


Lo = dxx + Vq ( x ) + no 


(4.8a) 

(4.8b) 
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At order e 3 we have 

L 0 u 3 = -AzzUg - AV 2 u e - a\A\ 2 A\u e \ 2 u e . 

The solvability condition of this equation is 

A zz -aA + a 1 \A\ 2 A = 0 1 (4.9) 

where 




V 2 uzdx, ct i = 


D 


\u e \ 2 u 2 e dx. 


(4.10) 

Equation m for the wave envelope A(Z) is our re¬ 
duced model for nonlinear wave dynamics near an ex¬ 
ceptional (phase-transition) point. Since V 2 and u e are 
■PT-symmetric and D real, a and ay are real. 

The reduced model (14.91) is a fourth-order dynamical 
system since A is complex. However it has two conserved 
quantities, 


h = \Az\ 2 ~a\A \ 2 + ^-\A\\ 


and 


I 2 = A*Az - AA 


Z ) 


(4.11) 


(4.12) 


where dlk/dZ = 0 (fc = 1,2). Due to these two conserved 
quantities, solution dynamics in Eq. (l4~9l) is confined to 
a two-dimensional surface, thus this dynamics cannot be 
chaotic. When Z — > oo, the solution A can only approach 
a fixed point, or a periodic orbit, or infinity (if oy > 0, 
infinity is further forbidden due to conservation of I\). 

The parameter a plays an important role in Eq. (14.91) . 
Let us consider the small-amplitude limit (|A| <C 1), in 
which case Eq. (HUD reduces to Azz — olZ = 0. If a < 0, 
these infinitesimal (linear) modes are bounded, meaning 
that the system is below phase transition. But, if a > 0, 
these linear modes exponentially grow, indicating that 
the system is above phase transition. Recall that a is de¬ 
pendent on the potential perturbation V 2 . Thus whether 
the perturbed potential is above or below phase transi¬ 
tion depends on the sign of a. In addition, the value of 
a also determines whether or not the underlying excep¬ 
tional point /r = /to is a phase-transition point: if a ^ 0, 
then n = n o is a phase-transition point; if a = 0, then the 
answer is not certain, and further analysis is needed in 
order to determine whether /r = /io is a phase-transition 
point or not. 

In the next two sections, we will describe the predic¬ 
tions of the reduced model (14.91) and compare them with 
the full system m- In all our numerical comparisons, 
we will use the Scarff-II potential (12.31) with Vr = 2. At 
this Vr value, an exceptional point occurs when 


Wo = 2.25, n o « —0.3144, 


(4.13) 


and this exceptional point is a phase-transition point. In 
the format (14.11) of the perturbed potential, this Scarff-II 
potential has 

Vq(x) = 2sech 2 (x) + i2.25sech(x)tanh(x), (4.14a) 

V 2 (x) = icsech(a:)tanh(a:). (4.14b) 


This potential is above phase transition when c = 1 and 
below phase transition when c = — 1. At this phase- 
transition point, the coefficients in the reduced model 
(USD are found to be 


a « 0.3144c, (7i « 0.2700cr. 


(4.15) 


For these coefficients the eigenfunction u e (x) has been 
normalized to have unit amplitude. In all our compar¬ 
isons, we take e = 0.2. This e is not very small, but 
predictions of the reduced model (14.91) still match those 
in the full system (12.11) as we will see below. 


V. SOLUTION BEHAVIORS ABOVE PHASE 
TRANSITION 

Our main interest is to investigate nonlinear wave dy¬ 
namics above phase transition (a > 0). Previous studies 
on nonlinear 'P'T-symmetric systems overwhelmingly fo¬ 
cused on solution behaviors below phase transition, be¬ 
cause it was argued that coherent structures such as soli- 
tons would be unstable above phase transition (at least 
in 'P'T-symmetric periodic potentials). We will show in 
this section that in PT-symmetric localized potentials, 
stable solitons and robust oscillating nonlinear modes do 
exist above phase transition. 


A. Soliton families and their stability 

First we consider soliton solutions, which correspond to 
constant-amplitude solutions in the reduced model (14.91) . 
Specifically, constant-amplitude solutions of the form 

A(Z) = A 0 e~ illlZ (5.1) 

in Eq. (14.91) correspond to soliton solutions of the form 


i/j = u{ x )e- i{ll0+ ^ 1 ' )z 


(5.2) 


in Eq. m, where u(x) ~ eAou e (x) to leading order. In 
this A-formula, Aq will be made real positive from phase 
invariance. Substituting (EH) into PI , we find /ii as 




(5.3) 


where the quantity under the square root must be non¬ 
negative. This equation relates the propagation constant 
/ii to the soliton amplitude parameter Aq. 

Since a > 0 above phase transition, solutions ex¬ 
ist only when > 0. For the Scarff-II potential (14.141) . 
this means that above phase transition, solitons can only 
exist under self-focusing nonlinearity (a > 0). The phys¬ 
ical reason for the existence of these solitons comes from 
the nonlinear feedback. It is commonly known that a 
■PT-symmetric complex potential is above phase transi¬ 
tion when the imaginary part of the potential (relative 
to the real part) is above a certain threshold. In the 





6 


current case, even though the linear potential V(x;e) is 
above phase transition, the nonlinearity-induced positive 
refractive index <r|'i/>| 2 , when added to this linear poten¬ 
tial, enhances its real part, which makes its imaginary 
part relatively weaker. As a consequence, the nonlinear¬ 
ity transforms the effective potential from above phase 
transition to below phase transition El- 

Notice also from Eq. (15.31) that these solitons exist only 
above a certain amplitude (or power) threshold, which is 


as 


Ag > a/a\. 


(5.4) 


This means that the nonlinearity-induced positive re¬ 
fractive index must be strong enough in order to trans¬ 
form the effective potential from above to below phase 
transition. Consequently, these solitons do not bifurcate 
from linear modes of the potential. In addition, the two 
branches of these solutions [corresponding to the plus 
and minus signs in (15.31) ] are connected at this amplitude 
threshold and thus belong to a single soliton family. 
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FIG. 2: (a) Power curve for the family of solitons above phase 
transition in the Scarff-II potential (14.141) with a = 1, c = 1 
and e = 0.2; solid blue indicates stable solitons and dashed 
red unstable ones; the horizontal dashed line is the analytical 
prediction for the power minimum; (b) profile of an exam¬ 
ple soliton at the black-dot point of the power curve (where 
= —0.5); (c) discrete eigenvalues for stable solitons versus 
the propagation constant fi; solid blue are numerical values 
and dashed red analytical predictions; the shaded region is 
the continuous spectrum; (d) the numerically obtained linear- 
stability spectrum for the soliton in panel (b). 


For the Scarff-II potential (14.14|) with a = 1, c = 1 and 
e = 0.2, we have numerically obtained these predicted 
solitons above phase transition, whose power curve is 
plotted in Fig.[2Ka). Here the soliton’s power is defined as 
P(p) = f_ | i/)\ 2 dx. It is seen that this numerical power 
curve indeed has a minimum threshold. The analytical 
prediction for this power threshold, obtained from Eq. 
(15.41) and the leading-order perturbation solution (14.71) 


P — 

1 m.rn. 


2 dx, 


CTl 


is also depicted in Fig.[2Ka) (as a horizontal dashed line). 
It is seen that this analytical power threshold matches 
the numerical value very well. At the black-dot point of 
the numerical power curve (where /i = —0.5), the pro¬ 
file of the corresponding soliton solution is illustrated in 
Fig. [21[b). This soliton is 'PT-symmetric, as are all other 
solitons in this family. 

Stability of these solitons can be analyzed by exam¬ 
ining the stability of constant-amplitude solutions E3) 
in the reduced ODE model (14.91) . Let us perturb this 
constant-amplitude solution by normal modes as 


A(Z) = (a 0 + Ae XAZ + P*e A ^ z ) 


,-ijUi Z 


where A, B <C 1, and A a is the eigenvalue from the en¬ 
velope equation. Plugging this into (14.91) and linearizing, 
we obtain 


La 


= 0, 


where 


a 


t _ ( \ 2 a — 2L\aPi + Pi T ct 

A f-if + a X 2 a + 2L\aPi + n 2 + a 


Pi 


Requiring the determinant of this matrix L A to vanish, 
non-zero eigenvalues Aa are then derived as 


A a = ±iy/ 2(3p 2 + a). 


Since a > 0 above phase transition, this formula predicts 
a pair of purely imaginary discrete eigenvalues, indicating 
that the constant-amplitude solution (15.11) is stable in the 
ODE model (14.91) . This implies that the soliton solution 
(15.21) is also stable in the original model (12.11) . Taking 
into account the scaling Z = ez, an approximation for 
non-zero discrete eigenvalues of this soliton is 


A ~ e\ A = 2(3p 2 + a). 


(5.5) 


Numerically we have confirmed the stability of these 
solitons near phase transition. This is achieved by com¬ 
puting the linear-stability spectrum of these solitons by 
the Fourier-collocation method [32} • For example, for the 
soliton shown in Fig. [2Kb), its linear-stability spectrum 
is displayed in Fig. [2Kd). All eigenvalues in this spec¬ 
trum are purely imaginary, indicating that the soliton is 
linearly stable. In addition, the pair of discrete imag¬ 
inary eigenvalues in this spectrum correspond to those 
predicted analytically by formula (15.51) . Similar compu¬ 
tations are performed for other solitons, and their stabil¬ 
ity is indicated by solid blue lines on the power diagram 
of Fig.[2Ka). Quantitative comparison between numerical 
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discrete imaginary eigenvalues and their analytical pre¬ 
diction (15.51) is made in Fig. H c )j and reasonable agree¬ 
ment can be seen (even though e = 0.2 is not small here). 

At high powers, we find that these solitons above phase 
transition become linearly unstable, and this instability 
is shown on the power curve of Fig. Ha) by dashed red 
lines. The instability on the left side of the power curve 
is induced by complex-eigenvalue bifurcations from edges 
of the continuous spectrum, while instability on the right 
side of the power curve is caused by complex-eigenvalue 
bifurcations from interiors of the continuous spectrum. 
These higlr-power solitons correspond to amplitude val¬ 
ues on the order A ~ 1/e and lie outside the validity of 
our perturbation theory, thus their instability does not 
contradict our stability result for low-power solitons. 


B. Oscillating solutions 

The behavior of solutions away from the soliton equi¬ 
libriums can be largely captured by focusing on the case 
where A is purely real in the reduced system (14.91) . In 
this case, the model equation becomes a simple second or¬ 
der ODE which we analyze using phase portraits. Above 
phase transition, a > 0, this breaks into two cases de¬ 
pending on the sign of the nonlinearity. 


1. Positive (Ji 

In this case, the phase portrait is shown in Fig. H a )> 
where a and ay values are taken from Eq. (14.151) with c = 
1 and <7 = 1 (focusing nonlinearity). This phase portrait 
contains three fixed points. One of them is the origin, 
which is unstable, signifying that the system is above 
pha se tra nsition. The other two fixed points are at A = 
± \Ja /ay, which are stable, and they correspond to the 
soliton of minimum power (with /i \ = 0) in Eqs. (15.31) - 
m- Away from these fixed points, the phase portrait 
features two types of periodic orbits which are separated 
by a figure-eight trajectory joined at the origin. Inner 
periodic orbits surround the non-zero fixed points, while 
outer periodic orbits undergo wider amplitude swings. 

These periodic orbits in the phase plane imply the ex¬ 
istence of robust oscillating solutions away from solitons 
in the full PDE (12.11) . and such oscillating solutions are 
confirmed in our direct evolution simulations of that sys¬ 
tem. To illustrate, two examples of such PDE solutions 
are displayed in Fig. Hb,c). Oscillations in panel (b) are 
stronger, and they correspond to outer periodic orbits 
in the phase portrait (a). Oscillations in panel (c) are 
weaker, and they correspond to inner periodic orbits in 
the phase portrait. This solution correspondence can be 
made more explicit by projecting the PDE solution onto 
the phase plane. To do so, we recall the perturbation 
solution (14.7I) - (14.8|I . which to order e 2 gives 

ip(x,z) = [eA(Z)u e (x) — i e 2 A'(Z)u g (x)] e~ , ' J ' oZ . (5.6) 



FIG. 3: (a) Phase portrait of the reduced model (1X91) with 
<7i > 0 above phase transition. (b,c) Full PDE simulations in 
the Scarff-II potential (14.141) with a = 1, c = 1 and e = 0.2 
under different initial conditions. In (a), solid blue lines are 
full PDE solutions in (b,c) projected onto the phase plane 
[with the outer curve for (b) and inner curve for (c)]. 


Taking the inner product of this equation with u g (x) and 
retaining only the leading-order term, we get 

1 r°° 

A(Z) = — / ^(x,Z/e)u g (x)dxe i ^ € . (5.7) 

J _ oq 

Taking the inner product of (15.61) with u e (x) and recalling 
the relation (14.51) . we get 

i r°° 

A\Z) = / ip(x, Z/e)u e (x)dx e IMoZ/e . (5.8) 

z-D J -oo 

In this way, the full PDE solution can be embedded in 
the phase portrait of the ODE model for comparison. As 
a technical matter, the projected quantities (A, A') from 
the PDE solution by (15.7I) - (15.8I) are complex in general. 
But we have found that if the initial condition of the PDE 
solution is chosen according to Eq. (15^1) . then the imagi¬ 
nary parts of the projected (A, A') remain very small for 
very long distances. Thus we neglect those small imagi¬ 
nary parts and plot only the real parts of the projected 
(A, A') in the phase plane. 

For the two PDE solutions in Fig. Hb,c), their phase- 
plane projections are displayed as solid blue lines in panel 
(a). It is seen that these PDE projections closely nrinric 
the periodic orbits of the ODE model. 

It is noted that these predictions of periodically- 
oscillating solutions in the PDE system are valid on the 
distance scale of z ~ 1/e. Beyond this distance scale, the 
PDE dynamics generally starts to deviate from the ODE 
predictions. Our numerics shows that over very long dis¬ 
tances, these oscillations in the PDE solution gradually 
intensify and eventually break up, which is caused by res- 













onance of nonlinearity-induced higher harmonics of these 
oscillations with the continuous spectrum in our opinion. 


2. Negative <j\ 

When fji < 0, the phase portrait is shown in Fig. [4] 
(left panel), where a and o\ values are taken from Eq. 
(14.151) with c = 1 and a = — 1 (defocusing nonlinear¬ 
ity). In this case, except for the origin (an unstable fixed 
point), all trajectories escape to infinity. Similar solution 
behaviors are observed in the full PDE (12.11) . An exam¬ 
ple is shown in the right panel of Fig. 01 where the PDE 
solution is seen to first decrease, and then rise to high 
amplitudes. The projection of this PDE solution onto 
the phase plane is displayed as a solid blue line in the 
phase portrait. This projection closely follows the tra¬ 
jectory of the ODE model. After the solution amplitude 
has reached the order A ~ 1/e (beyond the validity of 
our perturbation theory), PDE solutions can eventually 
saturate in amplitude while continuing to shed radiation 
and grow in power. 


0.5 

n n 

A; 0 

-0.5 


FIG. 4: (Left) Phase portrait of the reduced model (14.911 with 
<7 1 < 0 above phase transition. (Right) Full PDE simulation 
in the Scarff-II potential (14.141) with a = —1, c = 1 and 
t = 0.2. The solid blue line in the left panel is this full PDE 
solution projected onto the phase plane. 



VI. SOLUTION BEHAVIORS BELOW PHASE 
TRANSITION 


predicts that constant-amplitude solutions exist when 
M > y/[a[, i.e., soliton solutions exist when |/l/ — I > 
ey/[a[. In addition, the amplitude Aq (and hence power) 
of these solitons can be arbitrary. Numerically we have 
confirmed this prediction in the Scarff-II potential (14.141) 
with a = 1, c = —1 and e = 0.2. The numerically ob¬ 
tained power curves of these solitons are displayed in 
Fig. Eta). Stability of these solitons can be analyzed 
in the framework of the reduced model m, and the 
eigenvalue formula ED shows that these solitons are 
linearly stable, which agrees with the numerical findings 
in Fig. EJa) for solitons at lower amplitudes (where the 
perturbation theory is valid). At higher amplitudes, the 
solitons do become unstable, similar to the case above 
phase transition in Fig. Ha). 



FIG. 5: (a,b) Numerically obtained power curves for the fami¬ 
lies of solitons below phase transition in the Scarff-II potential 
(14.141) with c = —1, e = 0.2, a = 1 in (a) and a = —1 in (b). 
(c) Analytical prediction for the power curve and linear sta¬ 
bility of solitons in (b). In all of (a,b,c), solid blue indicates 
stable solitons and dashed red unstable ones, (d) Comparison 
of numerically obtained (solid blue) and analytically predicted 
(dashed red) unstable eigenvalues for solitons in (b). 


In this section, we consider the predictions of our re¬ 
duced model for solution behaviors below phase transi¬ 
tion, and compare them with PDE solutions. 


A. Soliton solutions 

Below phase transition, a < 0, the ODE model (14.91) 
admits constant-amplitude solutions ED for both signs 
of the nonlinear coefficient ay, meaning that solitons exist 
under both focusing and defocusing nonlinearities. But 
behaviors of solitons for the two signs of ay are very dif¬ 
ferent. 

When ay > 0, formula (15.31) . when rewritten as 

Al = {n\ + a)/cn, ( 6 . 1 ) 


When fji < 0, formulae ED and ED predict that 
solitons only exist in the propagation-constant interval 
of |/i — /l«o| < £\/M" w ith a limited range of amplitude 
values |A 0 | < sjot/ai. The analytically predicted power 
curve from Eqs. (1X71) . (1X51) and ED is 


PM 


(n - /r 0 ) 2 + ae 2 
Ol 



( 6 . 2 ) 


which is plotted in Fig. ED C )- Numerically we have ob¬ 
tained these solitons, whose power curve is shown in 
Fig. EDb)- This numerical power curve closely resembles 
the analytical prediction in (c). In particular, the exis¬ 
tence of a power upper bound is confirmed. This close 
agreement between the perturbation theory and direct 
numerics is understandable, since these solitons have low 
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powers and are thus within the regime of validity of the 
perturbation theory. 

The physical reason for limited power ranges of these 
solitons is that, under defocusing nonlinearity, if this 
power is too large, the negative nonlinearity-induced re¬ 
fractive index would transform the effective potential 
from below phase transition to above phase transition, 
rendering stationary solitons impossible. 

Stability of these solitons with limited power ranges 
can be analyzed in the framework of the reduced model 
(14.91) . In this case, the eigenvalue formula (15.51) predicts 
that these solitons are linearly unstable when 

l/r - no\ < e\/M/3, (6.3) 

and stable otherwise. This predicted instability and sta¬ 
bility is shown on the predicted power curve in Fig.[5jc). 
It is seen that solitons at the top part of the power curve 
are predicted as unstable and the bottom ones predicted 
as stable. Numerically we have determined the linear sta¬ 
bility of these solitons by computing their stability spec¬ 
tra, and the results are shown in Fig. EL). Clearly the 
numerical results match those of analytical predictions. 
Quantitatively we have also computed real eigenvalues of 
unstable solitons and plotted them in Fig. Ed), together 
with their analytical predictions in Eq. (ESI). This quan¬ 
titative comparison shows good agreement as well. 


B. Oscillating solutions 

Like the previous case above phase transition, robust 
oscillating solutions exist below phase transition as well. 
As before, we will unveil such solutions by focusing on 
the case of real A in the reduced model (14.91) . 

If cti >0, the phase portrait of the reduced model is 
shown in Fig. [S] (left panel), where a and ay values are 
taken from Eq. (14.151) with c = — 1 and a — i (focusing 
nonlinearity). In this phase portrait the origin is a stable 
fixed point, a reflection that the system is below phase 
transition. Surrounding the origin are periodic orbits of 
various sizes. This implies an abundance of robust oscil¬ 
lating solutions in the PDE system. Numerically we have 
confirmed the existence of these oscillating solutions, and 
an example is shown in Fig. [G] (right panel). Projection 
of this PDE solution onto the phase plane is plotted as a 
solid blue line in the left panel, and good agreement with 
the ODE orbit is seen. 

If ay <0, the phase portrait of the reduced model (14.91) 
is shown in Fig. Eta), where a and ay values are taken 
from Eq. (14.151) with c = — 1 and a = — 1 (defocusing 
nonlinearity). This phase portrait contains three fixed 
points: the origin which is stable, and A = ±\fa/<7\ 
which are unstable. The latter two fixed points corre¬ 
spond to the soliton with maximal power (at /i = /io) in 
Fig. Ec)- Away from these three equilibria, trajectories 
are divided into two categories: periodic orbits surround¬ 
ing the origin, and orbits which escape to infinity. Nu¬ 
merically we have found both types of solutions in the 
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FIG. 6: (Left) Phase portrait of the reduced model (14.91) with 
cti > 0 below phase transition. (Right) Full PDE simulation 
in the Scarff-II potential (14.141) with a — 1, c = — 1 and 
e = 0.2. The solid blue line in the left panel is this full PDE 
solution projected onto the phase plane. 



A 



FIG. 7: (a) Phase portrait of the reduced model (1X91) with 
<ti < 0 below phase transition. (b,c) Full PDE simulations 
in the Scarff-II potential (14.141) with rr — — 1, c = — I and 
e = 0.2 under different initial conditions. In (a), solid blue 
lines are full PDE solutions in (b,c) projected onto the phase 
plane [with the upper curve for (b) and lower curve for (c)]. 


PDE system (12.11) under the Scarff-II potential (14.141) 
with a = —1, c = —1 and e = 0.2, and two examples 
are displayed in Fig. Etb,c). Projections of the PDE so¬ 
lutions onto the phase plane in panel (a) indicate that 
the ODE model accurately describes the PDE dynamics. 


VII. SUMMARY AND DISCUSSION 

In this article, nonlinear wave propagation in VT- 
symmetric localized potentials was investigated analyt¬ 
ically near phase transition. Necessary conditions for 
phase transition were first derived based on a generaliza¬ 
tion of the Krein signature. Then rich nonlinear dynam¬ 
ics near phase transition was revealed through a multi¬ 
scale perturbation analysis, which yielded a nonlinear 
ODE model for the amplitude of the solutions. Above 
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phase transition, this ODE model predicted a family of 
stable solitons not bifurcating from linear modes under a 
certain sign of nonlinearity. In addition, it predicted per¬ 
sistent periodically-oscillating solutions away from soli¬ 
tons. Under the opposite sign of nonlinearity, it predicted 
unbounded growth of solutions. Below phase transition, 
solution dynamics was predicted as well. We have com¬ 
pared all analytical predictions with direct numerical cal¬ 
culations of the full PDE system and good agreement was 
obtained. 

The analytical results obtained in this article are help¬ 
ful for several reasons. First, it is known that phase tran¬ 
sition is a distinct and important phenomenon in PT- 
symmetric systems. Thus the analytical condition for 
phase transition in terms of PT-Krein signatures helps 
understand when phase transition can or cannot occur. 
Second, the analytical predictions of nonlinear dynam¬ 


ics near phase transition contribute to a global under¬ 
standing of solution behaviors in PT-symmetric systems. 
Thirdly, even though our analysis was performed only for 
the potential NLS equation m, a similar treatment can 
obviously be extended to other PT-symmetric systems 
near phase transition, and similar solution dynamics is 
expected in all such systems. 
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